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We present first results for a study with two mass-degenerate dynamical Chirally Improved (CI) 
fermions on lattices of spatial extent 2.4 fm. The CI Dirac operator obeys the Ginsparg-Wilson 
condition in good approximation. The pion mass we use is still large ff(A!Q MeV) for the 16 3 x 32 
lattices with lattice spacing of 0.15 fm. The hadron masses are obtained with the variational tech- 
nique and the results compared with earlier quenched calculations with similar lattice parameters. 
We find indications for the isovector, scalar meson ao(980) near the experimental value, in con- 
tradistinction to quenched results which always produced a mass value compatible with the first 
excitation ao(1450). 
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1. Introduction 

It will take some time before we will be able to compute hadron properties in full QCD simu- 
lations with dynamical quarks that have full Ginsparg-Wilson (GW) chiral symmetry. On the way 
to this aim we are working with the so-called Chirally Improved (CI) fermions. These are realized 
as truncated solutions to the GW equations for a general ansatz for the Dirac operator [jj], |2|]. For 
each site this fermion action includes several hundred neighbors with distances ranging up to three 
links. 

Extensive quenched calculations have demonstrated good chiral as well as good scaling be- 
havior (pfl. Within the Bern-Graz-Regensburg (BGR) collaboration we have obtained results for 
the hadron spectrum and mesonic low energy constants in quenched simulations. These included 
several lattice spacings and volumes and the three light valence quarks. Emphasis in these studies 
has been put on deriving sophisticated techniques to analyze excited hadron states [Q ||, ^] and we 
expect to utilize that experience for full QCD configurations eventually. 

The CI Dirac operator involves many terms and as such its implementation in an HMC program 
is technically non-trivial. First results were obtained on 12 3 x 24 lattices and presented in Ref. 

Hp. However, the small volume (linear spatial size 1.8 fm) did not allow reliable results for the 
hadron masses. 

Here we present first results for our HMC simulation on larger 16 3 x 32 lattices at lattice 
spacing a = 0.15 fm and for two different values of the quark masses. We use two mass degenerate 
light quarks, the Luscher-Weisz gauge action, and stout smearing. One HMC-trajectory has <^(100) 
steps for one unit of HMC-time. More technical details may be found in 

2. Equilibration 

The approach to equilibration has been monitored by standard tests, like tracing the configu- 
ration mean plaquette or the number of conjugate gradient steps in the Monte-Carlo accept/reject 
step. For these two quantities we find autocorrelation lengths around 2. We skip the first 100 con- 
figurations and then analyze every fifth. Table [I] specifies the parameters for the runs discussed 
here. 

On lattices of this size the determination of the pion mass is possible even for individual 
configurations. Figure [j] demonstrates the (satisfactory) time dependence of M n for every 5th con- 
figuration of, e.g., run B. 

run /3 LW cp n conf rc mea s a [fm] m AW i[MeV] M K [MeV] 
A 4.65 -0.06 425 65 0.153(1) 33(1) 473(4) 

B 4.70 -0.05 350 50 0.154(1) 41(1) 507(5) 



Table 1: Parameters of the two runs discussed: run sequence, gauge coupling, bare mass parameter co, 
number of configurations ra con f, number of configurations analysed n me as> lattice spacing a assuming that the 
Sommer parameter ro = 0.48 fm, AWI-mass, pion mass. 
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Figure 1: Time history of the pion mass determined from single configurations and pseudoscalar operator 
P = dy^u for run B. The quark smearing types used (narrow for tin and wide for ww) are indicated and show 
excellent agreement. The fits are cosh-fits to the plateau-range Af = 3 — 15. 



30 
25 
20 
15 
10 
5f 



0.04 0.05 

Re(A) 



Figure 2: Histogram of the smallest values of Re (A) (left) and real eigenvalues (right) for run A (in lattice 
units). The AWTmass for this run is at amp^l = 0.0264. 



We also determined the low lying eigenvalues of the CI Dirac operator (cf., || ^]). Figure || 
shows a histogram of the distribution of the smallest Re(A) and of the smallest real eigenvalues. 
They extend inside the GW unit circle, since the CI operator is not exactly chiral. 

The topological charge v is determined from the number of real modes counted according to 
their chirality. Tunneling is a well-known problem for the GW-exact overlap fermions and sophis- 
ticated techniques have been developed to overcome this obstacle JTfj| , 11]. Due to the diagonal- 
ization method that we used (ARPACK) we may miss a few of the modes further inside the unit 
GW-circle. For the CI operator we find frequent tunneling between sectors of different topological 
charge (cf., figure|3[). 

In [^] the eigenvalue distributions and properties as well as results for a hybrid approach (over- 
lap operator and CI configurations) are discussed in more detail. 

The axial Ward identity (AWI) allows one to define the renormalized quark mass through the 
asymptotic behavior of the ratio 



Z A (d t A 4 (p = 0,t)X(0)] 



Z,„2m = 2m^ 



(2.1) 



Z P (P(p = 0,t)X(0)) 

where X is any interpolator coupling to the pion and Z A , Z P and Z„, denote the renormalization 
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Figure 3: History (left) and histogram (right) of the topological charge v during the run sequence A; 100 
configurations have been discarded and then every 5th is plotted. 

factors relating to the MS-scheme at a scale of 2 GeV. These have been calculated for the quenched 



case at several values of the lattice spacing and came out close to one []12|]. Determination of the 
values for the dynamical case is under progress. For now we work with the ratio 



(d t A 4 (p = 0,t)P(0)) 
(P(p = 0,t)P(0)) 



2m AW i 



(2.2) 



defining the so-called AWI-mass. The effects due to smearing of the quark sources is accounted 
for by comparing with point sink correlators. 

In figure |] we show the pion masses obtained for run A vs. the AWI-mass. In addition to the 
single point where m va i = m sea we also exhibit three partially quenched data points, i.e., determined 
for different valence quark masses for the same dynamical sea quark configurations. We find 
excellent extrapolation to the chiral limit, following the leading order (GMOR) behavior. This 
is no surprise since we know that the leading order behavior is also found in purely quenched 
results. 
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Figure 4: Pion mass squared vs. AWI-mass for run A. The data point marked by the filled square is the 
value where the valence quark and the sea quark masses agree, the other three points are partially quenched, 
i.e., the valence quark masses differ from the sea quark mass. 
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Figure 5: Results for the ground state nucleon compared to the quenched results on a lattice with similar 
parameters [Q]. 

3. Hadron masses 

For the hadron propagators we applied the techniques discussed in [Q, [5|]. We used Jacobi 
smeared quark sources and the same set of interpolating operators for the hadrons. The correlation 
matrix between the different operators at source and sink were then analyzed invoking the so-called 
variational method. The presently available statistics did not allow to identify excited states and we 
therefore discuss here only results for the ground states. 

In the quenched case [[|, ||] the scale has been set with a Sommer parameter value of ro = 0.5 
fm. In full QCD one usually uses somewhat smaller values. Alternative methods use f K , M p or 
to set the scale. In most of our plots we exhibit the masses in units of the measured value of the 
nucleon mass. Only for the nucleon (and for the values of the lattice spacing in the table) we have 
to rely on a physical scale and there we use r$ = 0.48 fm. 

Figure || compares the results for the dynamical runs with quenched results obtained on lat- 
tices with approximately the same physical size and lattice spacing In figure ^| we present a 
comparison for the vector meson mass and for the negative parity nucleons; for the latter we can 
identify the lowest two states. In all cases the dynamical results lie slightly below the quenched 
ones. 

As a surprise came the behavior of the iso vector, scalar meson oq. In quenched calculations this 
state was always seen too high, compatible with the first excitation <3o(1450). In the dynamical runs 
we now see the lowest mass compatible with the ground state ao(980) (see figure [7]). This channel 
is coupled to Tir\ and we will have to study its momentum, volume and quark mass dependence for 
clear identification as a single particle state. 

4. Outlook 

We are currently extending our statistics for run A and plan to also produce gauge configura- 
tions at significantly smaller pion masses, for similar lattice parameters. This should allow us to 
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Figure 6: Left: Comparison of quenched (circles) and dynamical (squares) results for the p -meson mass 
in units of the measured nucleon mass. The filled squares denote fully dynamical results and the open ones 
partially quenched results, both for runs A and B. Right: We can identify the two lowest lying negative parity 
nucleons which extrapolate to the experimental values. 
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Figure 7: Contrary to the quenched runs we see in the isovector, scalar channel a low lying state compatible 
with the ground state oq. 



find at least the first excitations for the hadron states discussed. In particular the Roper and the 
Delta state are of high interest. Also the ao signal ought to be verified by further analysis. 
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